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Abstract 

Today’s modem aircraft is based on air-breathing jet propulsion systems, which use moving fluids as 
substances to transform energy carried by the fluids into power. Throughout aero-vehicle evolution, 
improvements have been made to the engine efficiency and pollutants reduction. 

The major advantages associated with the addition of ITB are an increase in thermal efficiency and 
reduction in NO x emission. Lower temperature peak in the main combustor results in lower thermal NO x 
emission and lower amount of cooling air required. 

This study focuses on a parametric (on-design) cycle analysis of a dual-spool, separate-flow turbofan 
engine with an Interstage Turbine Burner (ITB). The ITB considered in this paper is a relatively new 
concept in modern jet engine propulsion. The ITB serves as a secondary combustor and is located 
between the high- and the low-pressure turbine, i.e., the transition duct. 

The objective of this study is to use design parameters, such as flight Mach number, compressor 
pressure ratio, fan pressure ratio, fan bypass ratio, and high-pressure turbine inlet temperature to obtain 
engine performance parameters, such as specific thmst and thmst specific fuel consumption. Results of 
this study can provide guidance in identifying the performance characteristics of various engine 
components, which can then be used to develop, analyze, integrate, and optimize the system performance 
of turbofan engines with an ITB. 

Visual Basic program, Microsoft® Excel macrocode, and Microsoft® Excel neuron code are used to 
facilitate Microsoft® Excel software to plot engine performance versus engine design parameters. This 
program computes and plots the data sequentially without forcing users to open other types of plotting 
programs. A user’s manual on how to use the program is also included in this report. 

Furthermore, this stand-alone program is written in conjunction with an off-design program which is 
an extension of this study. The computed result of a selected design-point engine will be exported to an 
engine reference data file that is required in off-design calculation. 


Nomenclature 


/ Fuel/air ratio 

Cp Specific heat at constant pressure 

P t Total pressure 

T t Total temperature 

V Absolute velocity 
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M Mach number 

a Sound speed 

g c Newton’s constant 

m Mass flow rate 

R Gas constant 

T Temperature, installed thrust 

e Polytropic efficiency 

F Force uninstalled thrust 

S Uninstalled trust specific fuel consumption 

hpn Low heating value of fuel 

W Power 

a Bypass ratio 

y Ratio of specific heats 

r| Efficiency 

71 Total pressure ratio 

t Total temperature ratio 

Total enthalpy ratio 
r\ m Mechanical Efficiency 

Subscripts 

itb Inter-stage Turbine Combustors 

b Main burner 

d Diffuser 

n Nozzle 

r Ram 

/ Fan 

c Properties between upstream and main burner 

b Properties between main burner exit and ITB 
t Properties between exit ITB and downstream, total 

cL Low pressure compressor 

cH High pressure compressor 

tL Low pressure turbine 

tH High pressure turbine 

nf Fan-Nozzle 


1. Introduction 

This program is designed to give a parametric cycle analysis of an ideal Turbo Fan air-breathing 
propulsion system. In most common air-breathing propulsion engines, the “heart” of a gas turbine is the 
gas generator. It consists of three major components namely, compressor, combustor, and turbine as 
shown in figure 1.1. 

The idea behind gas generator is to convert intake air mixed with fuel into high temperature and high 
pressure gas. Depending on the applications of the gas turbine, the energy provided is extracted and used 
for different applications (turbojet, turbofan, turbo-shaft, turboprop, and ramjet) through different 
mechanisms. A turbojet engine can be constructed by adding an inlet and a nozzle as shown in figure 1 .2. 
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Figure 1.1. — A Gas Generator Propulsion System 



Figure 1.2. — Schematic Diagram for Turbojet Engine 

The nozzle converts internal energy of the hot gas into kinetic energy or a thrust. The work extracted 
by the turbine is used to drive the compressor. In the case of turbofan, turboprop, and turbo-shaft, the 
work from the turbine is required to drive a shaft for turbo-shaft, a fan for turbofan, and a propeller for 
turboprop in addition to drive the compressor. The ramjet engine consists of an inlet, a combustor, and a 
nozzle at the exit. It does not require the compressor because the inlet already has an air-compressing 
mechanism such that intake air has sufficient kinetic energy to increase its pressure. 

The main objective of this analysis it to determine the relationships between engine performance 
(primarily thrust F, thrust specific fuel S) to design parameters (compressor pressure ratio, fan pressure 
ratio, bypass ratio, etc.), to design constraints (burner exit temperature, compressor exit pressure, etc.), 
and to flight environment (Mach number, ambient temperature, ambient pressure, etc.). 


2. Aircraft Engine Parameters 

2.1 Performance Parameters 

2.1.1 Thrust . — Thrust is the force used to sustain the flight (thrust = drag), accelerate flight (thrust > 
drag), decelerate flight (thrust < drag). Using figure 3.1 for the control volume, we can apply momentum 
balance to the control volume. Uninstalled thrust F of a jet engine (single inlet and single exhaust) is 
given by 


(m 0 +m fuel )V e -m 0 V 0 

^ Ve M3 )^e 

Sc 


( 1 . 1 ) 


where m 0 ,riif are mass flow rates of air and fuel respectively 
Vo ? V e are velocities at inlet and exit respectively 
P Q , p e are pressure at inlet and exit respectively 
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For ideal case, the hot gas is expanded to the ambient pressure which gives P e = Pq. Equation (1.1) 
becomes 


F = 


(m Q + m fuel )V e - riiQ V 0 


£c 


( 1 . 2 ) 


The installed thrust T is given by 


T - F D in i et D noz (1.3) 

where D in i et and D noz are the drag force from the inlet and the nozzle. 

2.1.2 Specific Fuel Consumptions . — It is the rate of fuel use by the propulsion system per unit of 
thrust produced. The installed specific fuel consumptions, TSFC , and uninstalled specific fuel 
consumptions, S , are given by 


m / 

TSFC = — l 
T 


(1.4) 


5 


F 


(1.5) 


2.1.3 Efficiencies of an Engine . — Some of the following parameters will be used widely in this 
program namely, thermal efficiency, propulsive efficiency, and overall efficiency. 

The thermal efficiency characterizes the net energy output extracted (shaft work) from the engine 
divided by the available thermal energy (fuel). 


fim - 


Km 

Qin 


( 1 . 6 ) 


where, 


fir// 

thermal efficiency 

Kut = 

net power out of engine 

Qin = 

rate of thermal energy released ( riifh PR ) 


The propulsive efficiency defines the ratio between the engine power output and the power being used 
to run the aircraft. 


fip = 


T ■ Vc 


W, 


out 


(1.7) 
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where, 


OP 

propulsive efficiency of engine 

T 

thrust of propulsion system 

Vo 

velocity of aircraft 

W 0 ut = 

net power out of engine 


An overall performance of a propulsion system is given by the combination between thermal and 
propulsive efficiencies. 


^ O - 1 ! p * 1 TH 


( 1 . 8 ) 


where, 

t\q is overall efficiency. 


2.2 Notations for Compressible Flow 

Some useful quantity notation for compressible flow will be used in this report namely, stagnation 
temperature, stagnation pressure, and Mach number. 

Stagnation temperature or total temperature T t is defined as the temperature obtained when steadily 
flowing fluid is brought to rest adiabatically without involvement of works. Applying first law of 
thermodynamic to calorically perfect gas gives: 


h e = hf =h t =h + 


V 


( 1 . 9 ) 


where, 

h t is the enthalpy at stagnation condition. 

With the assumptions constant specific heat coefficient, the above equation can be written as: 


T t =T + 


2c r 


= T 


1 + 


y-1 


M z 


( 1 . 10 ) 


where, 



T t 

= 

stagnation temperature 

T 

= 

static temperature 

M 

= 

Mach number 

Y 

= 

ratio of specific heats 


Stagnation pressure or total pressure P t is defined as pressure reached when a steady flowing fluid is 
brought to rest adiabatically and reversibly. Using the isentropic relation, total pressure is given by 
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f 


y 

Y-l 


( 1 . 11 ) 


Pt=P 


' l + ^-P-M 2 ^ 


Ratio of total temperatures x and ratio of total pressure n across a component are denoted by 
subscript: d for diffuser, cL for low-pressure compressor, cH for high-pressure compressor, b for main 
burner, ITB for inter-stage turbine burner, tH for low-pressure turbine, tH for high-pressure turbine, n for 
nozzle, nf for fan nozzle, and / for fan. For example: 


total pressure leaving diffuser 
total pressure entering diffuser 

(1.12) 

total temperature leaving diffuser 
total temperature entering diffuser 

(1.13) 


2.3 Exceptions 

For the free stream, ram, we define x r as a ratio of total temperature to static temperature and n r as a 
ratio of total pressure to static pressure. 


x r =^ = \ + y -d-Ml 


(1-14) 




1 + - 


- Afn 


(1.15) 


A ratio between total enthalpy of the burner exit and ambient enthalpy, denoted by x^, is defined such 
that it will be one of the input parameters: 


x X-b = 


^ X—itb 


ifp Pf )t> urnerex it 

(c t) 

' P '0, ambient 
1 C P T t) ITB exit 


M, 


(1.16) 


(1.17) 


P /O, ambient 


2.4 Component Performance 

In this analysis, it is acceptable to assume that the working fluid in the engine can be idealized as a 
perfect gas. Properties of an ideal gas strongly depend on the temperature. This cycle allows fluid 
properties variation across the engine which assumes constant fluid properties from the main burner 
entrance upstream (c pc , y c ), from ITB entrance to main burner exit (c pU y r ), and from ITB exit downstream 

(■ Cpitb ? Y itb)- 

2.4. 1 Inlet and Diffuser. — Pressure losses occur due to the friction within the inlet wall. The pressure 
total ratio, 7ij, is always less than 1. 
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In supersonic flight, the pressure losses cause shock waves which produce greater pressure losses. 
The inlet total pressure is defined as the product of the ram pressure ratio and the diffuser pressure ratio. 
Therefore the portion of the pressure losses due to the wall friction and the shock waves is defined by: 


71 d ~ n d max ^1 r 


(1.18) 


From the Military Specification 5008B [1], the following relation is obtained: 


T| r =1 

IV =1-0.075 (M 0 -l) 


M o + 935 


for Mq < 1 
for l<M n <5 


for Mo > 5 


2.4.2 Compressor and Turbine. — Compressor is measured through two types of efficiencies, namely, 
isentropic efficiency and polytropic efficiency. Isentropic efficiency is defined by 

_ ideal work of compression for given n c ^ ^ 

actual work of compression for given n c 

-r(y-i)/y _i 


Polytropic efficiency is defined as 

_ ideal work of compression for a differential pressure change ^ ^ 

c actual work of compression for a differential pressure change 

With the assumption of constant e c , we can obtain the relation between x c and n c \ 

x c=7tc (7-i)/(ye c ) (1.2: 

Going through similar procedure as the compressor, we obtain turbine isentropic efficiency, turbine 
polytropic efficiency and relation between x t and n t as follows: 

_ ideal turbine work for given n t 


actual turbine work for given n t 

1 -t, 


_ T y/(y-i)^ 


(1.26) 
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3 4 4.4 4.5 5 


8 9 


Figure 3.1. — Station numbering of a turbofan engine with ITB. 


3. TurboFan — Separate Exhaust Streams with ITB Cycle Analysis 

The station numbering for the turbofan cycle analysis with ITB is in accordance with APR 75 5 A [2] 
and is given in figure 3.1. The ITB is located between station 4.4 and 4.5. 

The air stream entering the turbofan engine will flow through the fan and the engine core separately. 
Cycle analysis is then applied to both the bypass stream and engine core stream separately as listed below. 
The analysis follows closely as described in Mattingly [3]. 


3.1 Assumptions 


The following assumptions are employed: 

1 . Perfect gas upstream of main burner with constant properties y c , R c , c pc 

2. Perfect gas between station four and five with constant properties y t , R t , c pt 

3. Perfect gas downstream of inter-stage burner with constant properties y^, c p j t b 

4. All components are adiabatic 

5. The efficiencies of the high-pressure compressor ( HPC ), low-pressure compressor (, LPC ), 
fan, high-pressure turbine (HPT), and low-pressure turbine (LPT) are described through the 
use of polytropic efficiencies, i.e., e C H, e c i 5 ej, e t ^ and e t i, respectively. 


3.2 Fan Stream 

Step 1 . — Uninstalled thrust of fan stream, Fy; is given by 

F f =^{V l9 -V 0 ) + {P l9 -P 0 )A l9 


( 3 . 1 ) 
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Rearranging equation (3.1) gives: 


Step 2 . — 


Step 3 . — 


where 


Step 4 . — 


where 


F f ao|>19 

mj g c ^«0 V\ 9 /a o y c , 



(3.2) 


(3.3) 


(3.4) 


(3.4a) 


(3.5) 


(3.5a) 


Step 1 . — Uninstalled thrust 


3.3 Engine Core Stream 


F c - — (m 9 V 9 -m 0 Vq)+{P 9 -Pq)A 9 

gc 

Rearranging equation (3.6) gives: 


fc_ 

m 0 


a ()(m9 V 9 P « 

Sc \ ™C V ^9 y 


(3.6) 


(3.7) 
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Step 2. 


where 


^. = 1 + / 


m 9 =m c +m b + m itb 
f = fb + fitb 


A-*- 

m n 


r _ m itb 

Jitb 

m r 


Step 3 . — Multiplying 


a 9 p 9 


m,. 


V P 9J 


by 


Mg Uq CIq 

P9 A 9 V 9 Yc R c T 0 


and rearranging gives: 


a 9 p 9 


m r . 


1-^ 

V P 9 J 


= (1 + /)^ 


R itb T 9 /Tq a 0 


R c V 9/ a 0 Sc 

Step 4 . — Uninstalled thrust for the engine core becomes: 

F, 


\ P 9 J 


Y c 


c _ 

m c g c 


(i + /)i-i / 0 + (i + /)^M 

a 0 R c Vg/aQ 


( p ^ 

1-%- 

i 

v ^9 y 

Yc_ 


Step 5 . — 


fv ^ 2 

v _ 9 

V a 0 


= Mq 


Yjtb R itb T 9 
Yc R c p o 


From the total pressure and Mach number relation, 


Mo = 


<N 

( p \ 

1 

V P 9 y 


-1 


v 

t 9 _ /To 


L t 9 / 


'0 


To T<9 



(r,v*-l)/y lV6 


where 


(3.8) 

(3.8a) 

(3.8b) 

(3.8c) 

(3.8d) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 
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(3.13a) 


r /9 


- n r n d n cL n cH n b n tH n itb n tL n n 


T t9 _ 

x r x d x c x b X t fl Xj tb x tL x n 

h 


V = ! 


=1 


l t4 


T b = — 

T 0 T r t cH t cL 

Step 6 . — Apply First Law of Thermodynamics and ideal gas relation to the main burner: 

C pc Tft ^ b ^1 b hpR—b ~ ™4 ^ pt ^t4 


Multiply the above equation with — 


TtO 

™c c pc T 0 ^ t2 


, rearrange gives: 


T c T r + 


= (, + hW _ b 


C pc T 0 


Solve for fy. 


fb = 


^ X-b ^r^cL^cH 
P\bhpR pc^o)~ 


i-b 


where 


T X-b ~ 


C p t Tt4 

C pc T 0 


T to ~ T a (for the case adiabatic) 


T c - x cH x cL 


Step 7 . — Apply First Law of Thermodynamics and ideal gas relation to the ITB: 

™4 c pt T t4A + ™itb ^itb h PR-itb = ™4.5 c pitb T t4.5 
Rearranging the above equation gives, 


(3.13b) 

(3.13c) 

(3.13d) 

(3.13e) 

(3.14) 

(3.15) 

(3.16) 

(3.16a) 

(3.16b) 

(3.16c) 

(3.17) 
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(3.18) 


where 


■fa =(! + /») M~ T )T" — 

’H itb hpR '[Cpc^o) x X' - itb 


_ C pitb ^4.5 

^ A/ itb rri 

^ pc ^0 


T t9 ~ T i5 (for the adiabatic case) 


T c - t cH x cL 

Step 8. — Apply first law to each individual compressors and turbines: 

Power balance for LPC : W cL = m c c pc (T t2 5 -T t i) 

Power balance for HPC: W cH = m c c pc (T t3 - T t2 5 ) 

Power balance for HPT : W tH = m 4 c pt r\ m _ HP (T l4 - T t4A ) 

Power balance for LPT: W tL = m 4 5 c pitb r) m _ LP (T t4 5 - T t5 ) 

Power balance for fan: Wf =ihf c pc (/} j 3 ~T t2 ) 

It is chosen that the HPT and HPC are connected by a single shaft; therefore ideal turbofan the 
relation is given by: 

WcH = Wtff 

Multiplying both sides of equation (3.24) with — gives: 

™c C pc Tt2.5 


Solve for x^: 


where 


bcH - 0 - (l + fb) x cH x b 0 m-HP “ x tH ) 
C pc 


x tH - 1 + 


1 ~ x cH 

(l + fb) x cH x b T 1 m-HP 

c pc 


x b = 


T_t 4 

T 0 x r x c 


(3.18a) 

(3.18b) 

(3.18c) 

(3.19) 

(3.20) 

(3.21) 

(3.22) 

(3.23) 

(3.24) 

(3.25) 

(3.26) 

(3.26a) 


NASA/TM— 2005-2 13658 


12 



LPT ' LPC , and the fan are connected by a single shaft; therefore ideal turbofan, the relation is given by: 


W cL +W f =W tL 


Multiplying both sides of equation (3.27) with — 


Solve for T t p: 


where 


™c c pc ^tl 


gives: 


( t cL -l)+a(ly _ l) _ _P}}b_ (i _j_ y )_1^5_ x\ m _Lp (l — TfL ) 
C pc *t 2 


T tL ~ 1 + ‘ 


+a(l-x/) 


(1 + /) 


c j)itb T t 4.5 

c pc Tfl 


■H m-LP 


l t4.5 

Ttl 


~ x cL x cH x b x tH T itb 


L t 4.5 


^tib rri 

*0 t r t c Tfo X tH 

Step 9 . — Total uninstalled thrust per unit mass flow rate intake is given by: 


F c F f 


F m c 


m 


f 


rho (1 + a) 

Step 10 . — Thrust specific fuel consumption (5) is given by: 

e _ fb fitb 


(l + a )T~ 


m 0 


(3.27) 


(3.28) 


(3.29) 


(3.29a) 

(3.29b) 


(3.30) 


(3.31) 
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4. Summary of Equations 


All equations below are executed and computed in subroutine itb() inside Module 1 of Microsoft 
Visual Basic. The flow chart of this subroutine can be found in figure B.5 in appendix B. These equations 
are solved sequentially, except for MSH fuel bum model, which is solved by using successive substitution 
method. 

Inputs. T/q, 7q, y c , Cp C , Yt, Cp Y//£, Cpi fo, hpp.fo hpp_iffo Tlj ma x? ^itb> ^ nfi e cH-> e cL ? 

1\itb> r\m-HP, hm-LA ^V^9> ^0^1 9 ? T t4, ^4.5, ^ c //, a 

Outputs : F I rn , S,/ 

Equations : 


R =^ c ' c 

^ VC 

Yc 


(4.1) 

ii 

^ ", 

^ 1 


(4.2) 

p _ Y zYZ? ^ „ 
^zYZ? — c pitb 

Yitb 


(4.3) 

& 

II 

o 


(4.4) 

F 0 - M 0 a 0 


(4.5) 

t,=1+ Yc_1 M 0 2 

r 2 ° 


(4.6) 

"i 

II 

H 

rs 

Ci 


(4.7) 

>! 

II 

•-t 

VI 

o 


r| r = 1 — 0.075(M 0 -l) 1 ' 35 for 

1<M 0 <5 

(4.8) 

800 

Or = 4 for 

Mq + 935 

M 0 >5 


71 cf — 71 </ maxOr 


(4.9) 

_ Cpt 7*4 
T X-b ~ n T 
^ pc ^0 


(4.10) 

_ C pith T t 4 5 

itb /O rjl 

^ pc ^0 


(4.11) 

1cH=K C H 


(4.12) 
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(4.13) 


(4.14) 

x c = x cH x cL 

(4.15) 

■J iL 

£ 

II 

-C) 

P 

(4.16) 

T\cH ~ , 

- 1 

(4.17) 

j ^ 

*> ^ 
S J" 

O b-» 
P 

II 

13 

$r 

(4.18) 

iff -1 
0/ = 1 
V 1 

(4.19) 

If CSH model then 


j- _ x X-b ~ x r x cL x cH 

r \b h PR/(cpc T o)-'tX-b 

(4.20a) 

Else 


T t3 = T O x r x cL x cH 

(4.20b) 

ENTHALPY 0, hn) 


Set initial value of fuel/air ratio at station 4, 


/■ _ x X-b ~ x r x cL x cH 

J4i / \ 

Mb 11 PR / \ c pc T 0)~ x X-b 

(4.20c) 

3 ENTHALPY [T tA f 4i h t4 ) 


y _ hf4 ~ ht3 
VbhpR -ht4 

(4.20d) 

If //, - j\i > 0.0001, then / 4( - = f b and go to 3; else continue. 


End if 


T -1 1 l ~ x cH 

x tH~ [+ c 

(l + fb)— x cH x b 0 m-HP 
C pc 

(4.21) 
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If CSH model then 


fitb ~^ + fb) 


pkHtb pk^bHH 
'H itb^PR pc^o)~ ^X-itb 


Else 


T t4.4 ~ T t4 x tH 

ENTHALPY (T t4A , f b ,h t4A ) 

Set initial value of fuel/air ratio at station 4.4, 


f4.5i -(l + fb) 


^itb^PR /( C pcf) )~A-// 7 > 


4 ENTHALPY {T tA5 ,A 5i ,h t4 . 5 ) 


fitb ~{^ + fb) 


^4.5 ~ ^4.4 
Pitb h PR ~ h t4.5 


If | fb + fitb - >4.5*1 > 0-0001, then f 45i = f b +f itb and go to 4; else continue. 
End if 


(4.22a) 


(4.22b) 


(4.22c) 


(4.22d) 


f = fb + fitb 

(4.23) 

P 4 5 

Htb = T 

7q T r T c Tfo T t fj 

(4.24) 

T t 4 5 

* = x cL X cH H HH x ith 

HI 

(4.25) 

l-x ci +a(l-x f ) 

T t r = 1 + 

n i n Cpitb Tt4 - 5 n 

U + J ) T Pm-LP 

C pc *t2 

(4.26) 

%, =t„A 

(4.27) 

) e tL ] 

(4.28) 

_ 1 ~Hh 

^tH ~ ]/p 

1 -HH ,etR 

(4.29) 

n - l ~ TtL 

,L i W'*- 

(4.30) 

r,9 

— - n r n d n cL n cH n b n tH n itb n tL n n 
‘9 

(4.31) 
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(4.32) 



F J_ «o|>19 M() , T l9 /T 0 1-Pq/V 

™f V l9/ a 0 Y c y 


F_ 

m 0 



+ a 



m c nij- 


(1 + a) 


(4.33) 

(4.34) 

(4.35) 

(4.36) 

(4.37) 

(4.38) 

(4.39) 

(4.40) 

(4.41) 

(4.42) 

(4.43) 
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(4.44) 


g _ fb fitb 

(1 + a) — 
m 0 


Kut _ 

- — {l + /) • vi -V 0 2 +a- (r& - Vq )} 

(4.45) 

m c 

^ ’ 6C 


Qm = 

: fb ' h PR-b ' 4/) + fitb ' h PR-ITB ' 4 ITB 

(4.46) 

m c 

Kut /™c 
477/ - A / • 
Qin l m c 

(4.47) 


_F/m 0 {l + a)V 0 

Wout/fhc 

(4.48) 


4 0 = 4 TH ' 4 P 

(4.49) 


5. User’s Manual 

The Microsoft® Excel program is written in combination between spreadsheet neuron cells, visual 
basic, and macrocode. These three combinations provide user-friendly software that compilation and 
preprocessing are no longer necessary. The input is to well-labeled Microsoft® Excel cells. Most of the 
values are pre-specified and an example case is displayed. It is good practice to save this case by another 
name before modifying the spread sheet in order to retrieve the starting case. User obtains result plots 
right a way just by clicking some simple buttons. 


5.1 Definitions 

First of all, it is necessary for users to understand some icons that appear in the Microsoft® Excel 
sheets, for example, the Input sheet as shown in figure 5.1. Table 5.1 shown below lists down all the icons 
along with their descriptions. 


Table 5. 

. — -Different kinds of tools anc 

indicators 

Name 

Appearance (for example) 

Descriptions 

Combo box 


ITB 

ON 


Users change options through 
selecting the items in the pull down 
menu list. 

Button 

Calculate 


Clicking any button will execute the 
macrocode associated with it. 

User input 




This is where users input the values 
for variables. 

Indicator 1 




The values here change with 
altitude. Do not modify. 

Indicator 2 


All input values are read and shown 
here. Do not modify. 
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5.2 Program Outline 


The program mainly comprises seven types of sheets, namely CoverPage , Instructions , Input , plot 
sheets, data sheets, Singlept and Station. Meanwhile, the main code structure is shown in figure B.l in 
appendix B. 

5.2.1 CoverPage. — The CoverPage sheet contains the information of authors of this program. Any 
questions regarding the program can be addressed to us through email or by phone. 

5.2.2 Instructions. — First time users are strongly recommended to read this sheet before running the 
program. Since there are always possibilities of getting error computations such as division by zero, 
square root of a negative quantity, or over floating and under floating number, the program is written such 
that it will not crash if those errors are encountered during the computation. It will instead tell the user 
where the computation encounters those errors. In this sheet, you will find details of how to run the 
program and how to fix a problem if something goes wrong. 

This sheet also explains several assumptions made in the equations so that the users aware of some 
cases in the equations that have been idealized to simplify the problems. 

5.2.3 Input. — This sheet is where most of the inputs are specified. The program will check input 
values in this sheet to make sure that all the inputs are specified. It will tell the users if there are inputs 
that have not been specified. 

There are four Combo Box in the Input sheet ( Combo Box is a list box that displays a list of values 
and lets the users select one of the value in the list), namely ITB , Unit , Fuel Burn Model and Choose a 
Plot , as shown in figure 5.1. You need first to specify the option in these combo boxes, such as ITB , Unit 
and Fuel Burn Model , before moving on to combo box Choose a Plot. 

Combo Box Unit lets you specify the unit system. Currently, the program can handle two unit 
systems, i.e., British gravitational units (English) and International system of units (SI). 

Combo Box ITB lets you turn ON and OFF the Inter-stage Turbine Burner ( [ITB ). This feature 
provides a flexibility to choose two types of engine and they are engine with ITB-ON and engine with 
ITB-OFF. With this feature, you will be able to see how the engine performance behaves with ITB-ON 
and with ITB-OFF. 

Combo Box Fuel Burn Model let user choose two models for the computation of the fuel-air-ratio 
inside the burner, i.e., Constant Specific Heat ( CSH) model and the Modified Specific Heat {MSH) model 
[3,4]. 
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Once all the input values and the options in combo box ITB , Unit , and Fuel Burn Model are 
specified, you can specify an option in combo box Choose a Plot. This combo box provides several 
options to choose from, and they are: 

1 . Vs Compressor Pressure Ratio 

2. Vs Flight Condition Mach Number 

3 . Vs Fan Pressure Ratio 

4. Vs ByPass Ratio 

5. Vs Main Burner Exit Temperature 

In each option, it performs multiple calculations at different values of only one corresponding design 
variable. Once you select one of them, you will be directed to the relevant plot sheet. 

5.2.4 plot sheets. — Please refer to section 5.3 Multi-point Calculations for detailed information. 

5.2.5 data sheets. — Please refer to section 5.3 Multi-point Calculations for detailed information. 

5.2.6 Singlept. — Please refer to section 5.4 Single-point Calculation for detailed information. 

5.2. 7 Station . — Please refer to section 5.5 Engine Station for detailed information. 

5.3 Multiple-Point Calculations 

The following discussion describes the instructions on how to run the program for each option 
available in combo box Choose a Plot in Input sheet. Each option shares the similar code structure as 
shown in figure B.2 in appendix B. 

5.3.1 “Vs Compressor Pressure Ratio”. — When this option is chosen, users will be directed to a plot 
sheet, namely ST VS PIc , as shown in figure 5.2. This sheet is used to plot engine performances (e.g., 
specific thrust or thrust specific fuel consumption ) and other properties (e.g., overall efficiency, nozzle jet 
velocity ratios, overall efficiency, fuel-air ratio, ITB inlet velocity and temperature rise across ITB ) versus 
compressor pressure ratio (n c ) with different bypass ratio. In order to execute the code properly, users are 
required to follow those circled number 1 through 3 sequentially as shown in figure 5.2. 

First of all, user needs to specify the values for three input parameters (i.e., n c i, np and Mq) as well 
as the lower and upper bound limit for 44 Compressor pressure ratio ” (7i c ) and the increment (Att c ) (see 
fig. 5.3). For instance, the lower and upper bound limit are 12 and 50, respectively, whereas the 
increment is 1 . 



Figure 5.2. — 44 Vs Compressor Pressure Ratio” plot sheet 
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What follows is to compute the values for F/P-compressor pressure ratio (u c h). By clicking the 
4 Generate ’ button, a tabular data for the compressor (including k c l, n c //, and n c ) will be listed as shown 
in figure 5.4. The cyan cell to the right of n c indicates the total number of iteration for 7i c . The smaller the 
value for Att c , the larger the total number of n c \ and, therefore, more calculation at different values of 7i c . 

Thirdly, users are free to choose the number of bypass ratio (a) by selecting options from the pull 
down menu as shown in figure 5.5. If it is more than 3 bypass ratios, users are required to manually input 
any number greater than 3 by selecting ‘User-Defined’ option. 

At this point, users are ready to run the code and generate the plots, simply by clicking 4 Calculate ’ 
button. All computed data will be tabulated in seven separate data sheets, namely, DataSTPIc , 

DataSFPIc , OverallPIc , NVCorePIc , NVFanPIc , VitbPIc , and dTrisePIc. 

Users are advised not to modify or change values in any of these data sheets. Clearing data as well as 
plots can be done by clicking ‘Clear data & plot’ button (see fig. 5.5). 

To return to Input sheet for other runs with different design parameters, users simply click on 4 Input 
sheet ’ button next to 4 Calculate ’ button. 



Figure 5.3. — Input parameters in 44 Vs Compressor Pressure Ratio” plot sheet 



Generate 




^cL 

^cH 


28 

1.65 

7.272727 

12 

1.65 

7.378788 

13 

1.65 

8.434848 

14 

1.65 

9.090909 

15 

1.65 

9.69697 

16 

1.65 

10.30303 

17 

1.65 

10.90909 

13 

1.65 

11.51515 

19 

1.65 

12.12121 

20 

1.65 

12.72727 

21 

1.65 

13.33333 

22 

1.65 

13.93939 

23 

1.65 

14.54545 

24 


Figure 5.4. — Tabular data for the compressor pressure ratios 


NASA/TM— 2005-213658 


21 



Clear data & plot 

INPUT sheet 

Calculate 



01. = 

2 

1 d 



1 A|gh3 

Qti = 

4 



3 Alpha 

Q^ = 

6 

User defined 


Figure 5.5. — Pull-down menu for 
fan bypass ratios (a) 


5.3.2 “Vs Flight Condition Mach Number”. — When users choose this option, the program hides the 
plot sheets and data sheets from the previous use and shows a associated plot sheet (i.e., STVSMo ) and 
seven associated data sheets (i.e., DataSTMo , DataSFMo , OverallMo , NVCoreMo , NVFanMo , VitbMo , 
and dTriseMo .) Users will then be directed from Input sheet to ST VS Mo sheet in a second. Similarly, 
this sheet is used to plot engine performances (e.g., specific thrust or thrust specific fuel consumption) 
and other properties (e.g., overall efficiency, nozzle jet velocity ratios, overall efficiency, fuel-air ratio, 
ITB inlet velocity and temperature rise across ITB) versus flight Mach number (Mq) with different 
bypass ratio. 

Similar to the instructions for compressor pressure ratio, user first needs to specify the values for 
three input parameters (i.e., n c i , Tip and n c ) as well as the range of Mq and the increment (A Mq). 

Secondly, users are free to choose the number of bypass ratio (a) by selecting options from the pull 
down menu. At this point, users are ready to run the code and generate the plots, by simply clicking 
4 Calculate ’ button. All computed data will be tabulated in all seven separate sheets just mentioned above. 

5.3.3 “Vs Fan Pressure Ratio”. — Similar to flight Mach number, this option opens a plot sheet (i.e., 
STVSPIf) and seven data sheets (i.e., DataSTPIf DataSFPIf OverallPIf NVCorePIf NVFanPIf 
VitbPIf and dTrisePIf). It plots engine performances (e.g., specific thrust or thrust specific fuel 
consumption) and other variables (e.g., overall efficiency, nozzle jet velocity ratios, overall efficiency, 

ITB inlet velocity and temperature rise across ITB) versus fan pressure ratio (nfi with different bypass 
ratio. The instructions to run the program in this sheet are similar to those for ‘flight Condition Mach 
number’. (See section 5.3.2.) 

5.3.4 “Vs ByPass Ratio”. — This option opens a plot sheet (i.e., ST VS Alp) and seven data sheets 
(i.e., DataSTAlp , DataSFAlp , OverallAlp , NVCoreAlp , NVFanAlp , VitbAlp , and dTriseAlp). It plots 
engine performances (e.g., specific thrust or thrust specific fuel consumption) and other variables (e.g., 
overall efficiency, nozzle jet velocity ratios, overall efficiency, ITB inlet velocity and temperature rise 
across ITB) versus different bypass ratio (a) with different fan total pressure ratios. The instructions to 
run the program in this sheet are analogous to those described above. The only difference is that selection 
of a is replaced by np 

5.3.5 “Vs Main Burner Exit Temperature”. — This option opens a plot sheet (i.e., ST_VS_Tt4) and 
seven data sheets (i.e., DataSTt4 , DataSFTt4 , OverallTt4 , NVCoreTt4 , NVFanTt4 , VitbTt4 , and 
dTriseTt4). It plots engine performances (e.g., specific thrust or thrust specific fuel consumption) and 
other variables (e.g., overall efficiency, nozzle jet velocity ratios, overall efficiency, ITB inlet velocity and 
temperature rise across ITB) versus main burn exit temperature (ff) with different flight Mach number. 
The instructions to run the program in this sheet are similar to those for ‘flight Condition Mach number’. 
(See section 5.3.2.) 
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5.3.5 Summary. — In summary, you need to do the followings to run the program: 

1 . Go to the Input sheet 

2. Specify the Unit system (SI or English) 

3. Specify the ITB switch. Is the computation for engine with ITB-ON or for engine with ITB-OFF? 

4. Specify the Fuel Burn Model system (CSH or MSH) 

5. Enter all the input parameters indicated in green cells (do not modify or change the value 
indicated in cyan). 

6. Specify Choose a Plot. You will be directed to a new sheet depending on the selection. 

In sheet STVSPIc (Compressor Pressure Ratio), 

a. Specify number of bypass ratio (a). 

b. Specify all input parameters in green cells. 

c. Click 4 Generate ’ to generate high-pressure compressor total pressure ratio. 

d. Click 4 Calculate ’ to compute and plot the results. 

e. Repeat the above steps for different input parameters. To return to Input sheet, simply click 
4 Input sheet ’ button. 

In sheet ST_ VS_Mo_( Flight Condition Mach Number), 

a. Specify number of bypass ratio (a). 

b. Specify all input parameters in green cells. 

c. Click 4 Calculate ’. 

d. Repeat the above steps for different input parameters. To return to Input sheet, simply click 
4 INPUT sheet ’ button. 

In sheet STVSPIf ( Fan Pressure Ratio), 

a. Specify number of bypass ratio (a). 

b. Specify all input parameters in green cells. 

c. Click 4 Calculate ’ to compute and plot the results. 

d. Repeat the above steps for different input parameters. To return to Input sheet, simply click 
4 Input sheet ’ button. 

In sheet ST_ VS Alp (ByPass Ratio), 

a. Specify number of fan total pressure ratio (nj). 

b. Specify all input parameters in green cells. 

c. Click 4 Calculate ’. 

d. Repeat the above steps for different input parameters. To return to Input sheet, simply click 
4 INPUT sheet ’ button. 

In sheet ST_VS_Tt4 (Main Burner Exit Temperature), 

a. Specify number of flight Mach number (Mo). 

b. Specify all input parameters in green cells. 

c. Click 4 Calculate ’. 

d. Repeat the above steps for different input parameters. To return to Input sheet, simply click 
4 INPUT sheet ’ button. 
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5.4 Single Point Calculation 


All instructions described above are termed as the multi-point calculations at different values of only 
one design parameter (e.g., n c , Tip Mq , a, or T P 4) The computed results show how the performance of a 
family of engines was determined by engine design choices and flight Mach number. 

Once the most promising engine is selected by users (i.e., all engine components are fixed, e.g., tt c , Tip 
altitude, a etc.), single point calculation is then run to get the engine performance at this particular 
operating conditions. These output data represents the selected engine’s design point condition. In order 
to predict the engine performance at off-design conditions, these outputs are then used as the reference 
conditions in an off-design code. This step is important because different reference conditions (i.e., 
different engines) will give different performance over the whole flight spectrum. 

For user’s reference, the flow chart of this computation is shown in figure B.3 in appendix B. 

5.4.1 Exporting Reference Engine Data . — For convenience, authors include Tnput/Output files' 
feature in this code. The output data from this on-design code can be exported to be a reference engine 
data file, which is required in off-design calculation. As a result, users do not need to manually re-input 
the reference conditions in off-design code, which is a time-consuming process. 

Users are to follow the step 1 to 13 to generate a reference engine data file: 

1 . Select Input sheet (as shown in fig. 5 . 1 ). 

2. Input the engine characteristics in Input sheet. 

3. Select SinglePt sheet (as shown in fig. 5.6). 

4. Input values of independent variables (i.e., Mq, a, altitude, n c and throttle settings) in green cells. 

5. Click on 4 Perform calculation ’ button to calculate the selected engine performance at a design 
point. 

6. Both the input values from Input sheet and the solution will be shown in light-blue cells. 

7. Click on 4 Export reference engine data ’ button to export these reference engine data to a text 
file. 

8. Figure 5.7 will appear and request for specifying the path and filename. 

9. Users may need to enter a new path if necessary, for example, 46 c:\itb\reference.txf\ The default 
path and filename is “c:\reference.txf \ Figure 5.8 appears when one of these situations happens: 

(i) The specified directory (default drive is 6 c:\’) does not exist. 

(ii) The file is being used by another program. 

10. Check for existing drive path or/and specify another filename. (Restart from step 3 again). 

1 1 . All reference data and engine characteristics of a specific design-point engine are successfully 
exported. 

12. Repeat step 1 to 1 1 for different engines. 
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On-design Single Point Calculation 


Inout check: 

Inout Desian 1 

variables: 





Pd^P-id _ 

1 

Mo = 

0 

h = 

0 

km 


19 = 

1 

a = 

8.4 

Tn = 

288.15 

K 




KcL = 

1.3 

Pn = 

101 325.00 

Pa 


mm = 

0.99 

Hc = 

39.3 

Perform calculation 

INPUT sheet 


0.95 

Hf= 

1.3 




J 

^ITH = 

0.95 

Ti+ = 

1600.0 K 


Export reference enqine data 

Hn = 

0.99 

T i+5 = 

1300.0 K 





■jw= 

0.98 

T k£-rTHoff = 

1300.0 K 







Outputs: 





e cH = 

0.9085 

K IH = 

0.1 52272 

11 IH = 

0.91 7937314 


Remarks: 

e cL = 

0.8738 

Hi L = 

0.467057 

11 IL = 

0.926882611 


■ => ln| 

e IH = 

0.8999 

11 CH = 

30.23077 

11cH = 

0.857801144 


=> Dc 

e, L = 

0.9204 

H C L = 

1.3 

11cL = 

0.8690081 39 


=> Dc 


0.93 

7I r = 

1.3 

flr = 

0.92733971 





r IH = 

0.68216 





rib = 

0.99 

$ IL = 

0.84606 

P|1^Pl9 = 

1.261 348375 



*1 rrB = 

0.99 

T CH = 

2.928886 






0.92 

r cL = 

1 .089799 

T,i^n = 

1.0841721 54 



11 m-LP = 

0.93r 1 

Tr 

1.08415 

Vis^n = 

58.95568433 



hpR-i = 

43124 kJ /kg 

fuel model = 

MSH 

v 9 /v n = 

210.5988768 



hpRvw = 

431 24 kJ/kg 

1 b = 

0.020891 

wr D = 

3.075658578 





* rre = 

0.006565 

PbPb = 

2.4724601 07 



Cpc = 

1.001994 kJ/kg-K 







c Pf = 

1.270719 kJ/kg-K 

M s = 

1 

m 5 = 

1.240438015 



^pit* = 

1 .20289 kJ/kg-K 

Mi B = 

0.585351 

M,b = 

0.585350632 



Tc = 

1.401343 

ffi = 

8.4 





7t = 

1.291711 

FAndat = 

254.2952 (N-sec/Kg) 




T» = 

1.31 3314 

S = 

11.48572 (mg/N-sec) 





rip = 1 .878652 % 
ti th = 36.95226 % 
rio = 0.694204 % 


Figure 5.6. — On-design single point calculation sheet. 



Figure 5.7. — Dialog box for specifying the path and filename 
where the reference conditions data will be stored. 
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Microsoft Excel 




There was a problem writing "h:\reference.txt". Please input another filename or close the application using this file, 


OK 


Figure 5.8. — This dialog box appears when there is an error 
writing reference conditions data to a specified file. 


Engine Station Data 


Inputs 


Mo = 

0 


CL = 

3.4 


n LPC - 

1.3 


■nc = 

39.3 


Hf = 

1.3 


Tt4 = 

1972.22222 

K 

Tt4.5 = 

1622.22222 

K 

h = 

0 

km 

t d = 

233.15 

K 

Pd - 

101325.00 

Pa 

Outnut 



ST = 

29.71 


TSFC = 

0.50 


nt = 

44.64 

% 

n P = 

1.27 

% 

no = 

0.57 

% 

f= 

0.033597 


4 = 

0.034203 


fitb — 

0.004394 



Station 

T. 

r 

T 

P 

M 

0 

513.63 

14.69 

513.67 

14.69 

0.00 

2 

513.63 

14.55 

493.93 

12.26 

0.50 

2.5 

565.14 

13.91 

547.62 

16.94 

0.40 

3 

1650.99 

571.67 

1621.79 

537.06 

0.30 

4 

3550.00 

543.03 

3114.04 

543.03 

1.00 

4.4 

2694.07 

133.72 

2694.07 

133.72 

0.30 

4.5 

2694.07 

127.03 

2353.05 

69.67 

1.00 

5 

2339.26 

70.55 

2307.07 

60.25 

0.50 

9 

2339.26 

69.34 

1690.36 

14.69 

1.70 

13 

562.22 

13.91 

544.79 

16.94 

0.40 

19 

562.22 

18.53 

526.13 

14.69 

0.59 


Update Station Data 

INPUT sheet I I 


Figure 5.9. — Engine station data sheet 


5.5 Engine Station Data 

The last sheet in the program is the Station sheet. This sheet displays the engine station flow 
properties (i.e., temperature and pressure) in forms of a table (see fig. 5.9) and charts (see fig. 5.10) based 
on the user input Mach number data at certain stations. 

Most of the input parameters are from the Input sheet. Additional inputs are also needed, such as 
those listed in “Inputs” column (i.e., to the left of the engine station data table) and Mach numbers at 
some engine stations (as shown in fig. 5.9). To click the 6 Update Station Data 5 button will update the 
flow properties at each engine station. Please refer to figure B.4 in appendix B for the code structure of 
this sheet. 
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3500 



(a) Temperature 



(b) Pressure 

Figure 5.10. — Examples of pressure and temperature vs engine stations 

5.6 Discussions 

A typical value for local Mach number at the ITB entrance is selected to be 0.3. This value 
( M 4 4 = 0.3) is then used to compute ITB inlet velocity (V44), on which one of the factors ITB design 
was based. In fact, each plot sheet will show the change of ITB inlet velocity with the corresponding 
design parameter, for example, V4 4 versus Mq in STVSMo sheet. 

At ITB-OFF condition, a plot of “Temperature Rise across ITB ” in each plot sheet (e.g., STVSPIc , 
ST VS Mo , STVSPIf, and ST VS Alp) will be cleared out. It is because the ITB-OFF engine acts as if 
it is a turbofan engine without any increase across ITB. 

5.7 Troubleshooting 

This program has been debugged several times. Therefore, whenever you encounter computation 
errors due to either zero division or square root of a negative quantity, you will be notified by a pop up 
window indicating where the computation problem is. It will tell you on which equation that the 
computational error is encountered. If you have any comments or bug problems you encounter in the 
program, you can report them to us for further improvement. Details about the contact number can be 
found in the CoverPage sheet of the program. 
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ITB Cycle Analysis: Warning / Error IX 


(a) 



ITB Cycle Analysis: Warning / Error 



(b) 

Figure 5.1 1. — Pop-up windows showing warning or error messages: 
(a) Equations where the error occurs; (b) The current values of 
design variables while error is encountered. 


5.7.7 Out of Range . — In multiple-point calculation, users may always see these pop-up windows as 
shown in figures 5.11a and 5.1 lb. There are two possibilities of having this error or warning messages: 


(i) The input values are incorrect. For instance, the main burner exit temperature (T t4 ) may be 

too small for an engine to produce enough thrust. Please double check all input values and try 
again. 

(ii) During the multiple-point computation, the design variable may be out of range. It will 
happen, especially for these design variables such as flight Mach number, fan pressure ratio 
and bypass ratio. For instance, the specific thrust of an engine is approaching zero as the 
flight Mach number is 3. If the upper limit of flight Mach number is set to be greater than 3 
(e.g., 5), the code will not be able to proceed at Mach number higher than 3 and thus give 
these warning messages. 

5. 7.2 Plotting errors. — It is also predicted that the plotting macro can encounter some problems in the 
future if the users do not fully understand the program. Therefore, authors have set up a way in order to 
fix the problem. In case the plot starts giving problem, you need to do the followings (this apply to all plot 
sheets: ST_VS_PIc , ST_VS_Mo , ST_VS_PIf ST_VS_Alp, and ST_VS_T t4 ): 

a. Go to the plot sheet (e.g., STVSMo) 

b. Go to cell D5 

c. Note that you will not be able to click on that cell because it is lying on the back of the plot chart. 

You may want to use arrow button from your keyboard in order to move to cell D5. 

d. In that cell, you will see “SUCCESS”. Change the character to “FAIL”. 

e. Clear the series in the chart plot by clicking right mouse button (RMB)on the chart plot. 

f. In RMB menu list, choose Source Data 

g. In ‘Source Data’ windows, pick Series tab and remove all the series. 

Users are not expected to understand macro code in the program. However the part that takes care the 
computation can be found in Modulel Sub Itb{ ). To open Microsoft® Excel Visual Basic windows, users 
need to press Alt + Fll. Any modification can be made in subroutine Itb{ ). Do not add new variables 
into the program! 
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Appendix A — Macro Visual Basic Code 


Sub itb() 

Dim rc, rb, rt, aO, vO, etar, taualpb, taualpitb 

Dim tau_c, taub, tauitb 

Dim T19T0, fe me, ffmf, denu, nume, vl9, v9 

Style = vbYes 
Style 1 = vbYes 

' Msg = "Do you want to jump to the next alpha?" 

Msgl = "Division by zero => Check zero inputs for gama's (eqn 1,2,3)" 

Msg2 = "Imaginary sound speed ==> Check negative inputs for gamac, gc or TO (eqn 4)" 

Msg3 = "Division by zero ==> gamac cannot = 1 (eqn 7)" 

Msg4 = "Division by zero => Check zero inputs for epe, TO, ecH, ecL & ef (eqn 10, 1 1, 12, 13 & 14)" 

Msg5 = "Division by zero ==> Check zero inputs for tau r & tau c (eqn 16)" 

Msg6 = "Division by zero => tau cH, tau cL & tau f cannot =1 (eqn 17, 18, 19)" 

Msg7 = "Invalid (negative) tau tH ==> Refer to equation 21" 

Msg8 = "Imaginary M9 ==> Pt9/P9 < 1 (eqn 31 & 34)" 

Msg9 = "Imaginary M19 => Ptl9/P19 < 1 (eqn 36 & 37)" 

Msg 10 = "Division by zero => denominator of 'fb' eqn cannot =0 (eqn 20)" 

Msgl 1 = "Division by zero => denominator of 'fitb' eqn cannot =0 (eqn 22)" 

wdata = 1 

fbi = 0 

fitbi = 0 

'warn = 0 

'If gamac - 0# Or gamat - 0# Or gamaitb - 0# Then 
'response = MsgBox(Msgl, style, Title) 

'Msgl = "Occurs at " & "PIc=" & pi_cH * picL & ", PIf=" & pi_f & ", Mo=" & M0 & ", alpha=" & alpha 
'response = MsgBox(Msgl, Style, Title) 

'If singlept <> 1 Then jump = MsgBox(Msgsheet, Style 1, Title) 

'GoTo 10 
'End If 

'added on 5/13/04 for controling ITB switch -1 
itbinit = 1 

If Sheets("Input").ComboBox2.Value = "OFF" Then itbinit = 0 
'added on 5/13/04 for controling ITB switch -2 
'Constants 

Call enthalpy(l, T o, 0, 0, Cpc, gamac) 

Call enthalpy(l, Tt4, 0, 0, Cpt, gamat) 
rc = (gamac - 1#) * Cpc / gamac '(1) 
rt = (gamat - 1#) * Cpt / gamat '(2) 

If gamac < 0# Or rc < 0# Or T o < 0# Or gc < 0# Then 
response = MsgBox(Msg2, Style, Title) 

'Msg2 = "Occurs at " & "PIc,PIf, Mo, ByPass = " & pi cH * picL & "," & pi f & "," & M0 & "," & alpha 
'Msg2 = "Occurs at " & "PIc=" & pi_cH * pi cL & ", PIf=" & pi f & ", Mo=" & M0 & ", alpha=" & alpha & 
", Tt4=" & Tt4 

'response = MsgBox(Msg2, Style, Title) 

'warn = 1 

'If singlept <> 1 Or last = 0 Then jump = MsgBox(Msgsheet, Stylel, Title) 

GoTo 10 
End If 

aO = (gamac * rc * T o * gc) A 0.5 '(4) 
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vO = aO * MO 


'( 5 ) 


If gamac =1# Then 

response = MsgBox(Msg3, Style, Title) 

’Msg3 = ’’Occurs at ” & ”PIc=” & pi_cH * picL & ”, PIf=” & pi_f & ”, Mo=” & MO & ”, alpha=” & alpha & 
”, Tt4=” & Tt4 

’response = MsgBox(Msg3, Style, Title) 

’warn = 1 

’If singlept <> 1 Then jump = MsgBox(Msgsheet, Style 1, Title) 

GoTo 10 
End If 

tau r = (gamac - 1#) * (0.5 * MO A 2) + 1# ’(6) 
pi_r - tau r A (gamac / (gamac - 1)) ’(7) 

If MO <= 1# Then ’(8) 

etar =1# 

End If 

If MO > 1 And MO <= 5 Then 
eta_r = 1# - 0.075 * (MO - 1#) A 1.35 
End If 

If MO > 5 Then 

eta r = 800# / (935# + MO A 4) 

End If 

pi_d - pidmax * eta_r ’(9) 

If picL < pi_f Then picL = pi_f 
picH - pi c / pi cL 

If Cpc - 0# Or T o - 0# Or gamac - 0# Or e_cH - 0# Or e_cL - 0# Or e_f = 0# Then 
response = MsgBox(Msg4, Style, Title) 

’Msg4 = ’’Occurs at ” & ”PIc=” & pi_cH * pi cL & ”, PIf=” & pi f & ”, Mo=” & MO & ”, alpha=” & alpha & 
”, Tt4=” & Tt4 

’response = MsgBox(Msg4, Style, Title) 

’If singlept <> 1 Then jump = MsgBox(Msgsheet, Style 1, Title) 

’warn = 1 
GoTo 10 
End If 
24 

tau alp b = Cpt * Tt4 / (Cpc * T o) ’(10) 

tau cH = picH A ((gamac - 1#) / (gamac * e_cH)) ’(12) 
tau cL = pi cL A ((gamac - 1#) / (gamac * e_cL)) ’(13) 
tau f = pi_f A ((gamac - 1#) / (gamac * e_f)) ’(14) 

tau_c - tau_cH * tau_cL ’(15) 

If tau_r - 0# Or tau c - 0# Then 
response = MsgBox(Msg5, Style, Title) 

’Msg5 = ’’Occurs at ” & ”PIc=” & pi_cH * pi cL & ”, PIf=” & pi f & ”, Mo=” & MO & ”, alpha=” & alpha & 
”, Tt4=” & Tt4 

’response = MsgBox(Msg5, Style, Title) 

'warn = 1 

’If singlept <> 1 Then jump = MsgBox(Msgsheet, Style 1, Title) 

GoTo 10 
End If 

tau b = Tt4 / (T o * tau r * tau c) ’(16) 
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If tau cH = 1# Or tau_cL = 1# Or tauf = 1# Then 
response = MsgBox(Msg6, Style, Title) 

'Msg6 = ’’Occurs at ” & ”PIc=” & pi_cH * picL & ”, PIf=” & pi_f & ”, Mo=” & MO & ”, alpha=” & alpha & 
”, Tt4=” & Tt4 

’response = MsgBox(Msg6, Style, Title) 

’warn = 1 

’If singlept <> 1 Then jump = MsgBox(Msgsheet, Style 1, Title) 

GoTo 10 
End If 

eta cH = ((taucH A e_cH) - 1#) / (tau cH - 1#) ’(17) 

eta cL = ((tau cL A e_cL) - 1#) / (taucL - 1#) ’(18) 

eta f = ((tau_f A e_f) - 1#) / (tau_f - 1#) ’(19) 

’add MSH on 7/30/04 -1 
'Component: MAIN BURNER 
H O = Cpc * T o 

tempfb - etab * hPRb / h O - taualpb 
If temp fb - 0# Then 
response = MsgBox(MsglO, Style, Title) 

’Msg7 = ’’Occurs at ” & ”PIc=” & pi_cH * picL & ”, PIf=” & pi f & ”, Mo=” & MO & ”, alpha=” & alpha & 
”, Tt4=” & Tt4 

'response = MsgBox(Msg7, Style, Title) 

'If singlept <> 1 Then jump = MsgBox(Msgsheet, Style 1, Title) 

'warn = 1 
GoTo 10 
End If 

’Fuel model 1:CSH 

If Sheets(”Input”).ComboBox4. Value = "CSH” Then 
fb - (tau alp b - tau r * tau c) / temp fb '(20a) 

Else 

'Fuel model 2: MSH 
loop3 - 0 

Tt3 = T o * tau_r * tau c ’eqn(20b) 

Call enthalpy(0, Tt3, 0, h_t3, 0, 0) 

f_4i = (tau alp b - tau_r * tau c) / temp fb ’eqn(20c), initial value for f_4i 

3 

If loop3 > 20 Then 
GoTo 10 
End If 

Call enthalpy(0, Tt4, f_4i, h_t4, 0, 0) 
temp fb = eta b * hPRb - h_t4 
If temp fb - 0# Then 
response = MsgBox(Msgl0, Style, Title) 

’Msg7 = "Occurs at " & "PIc=" & pi_cH * pi cL & ”, PIf=” & pi_f & ”, Mo=” & M0 & ”, alpha=” & alpha 
& ”, Tt4=" & Tt4 

'response = MsgBox(Msg7, Style, Title) 

'If singlept <> 1 Then jump = MsgBox(Msgjump, Stylel, Title) 

'warn = 1 
GoTo 10 
End If 

fb - (h_t4 - h_t3) / temp fb ’eqn(20d) 

If Abs(fb - f_4i) > 0.0001 Then 
f_4i = fb 

loop3 = loop3 + 1 
GoTo 3 
End If 
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End If 

’lkhl - added cp=cp(T,f) and gama=gama(T,f) on 1 1/15/04 
Errfb = Abs(fbi - fb) 

If Err fb > 0.0001 Then 
fbi = fb 

Call enthalpy(l, Tt4, fb, 0, Cpt, gamat) 

GoTo 24 
End If 
'lkh2 

’add MSH on 7/30/04 -2 

tautH = 1# - (taucH - 1#) * taucL * taur / ((1# + fb) * taualpb * effmtH) '5/1 1/04 

’(21) 

If tau tH < 0# Then 

response = MsgBox(Msg7, Style, Title) 

’Msg7 = ’’Occurs at " & "PIc=" & pi_cH * picL & ", PIf=” & pi_f & ", Mo=" & M0 & ", alpha=" & alpha 
& ", Tt4=" & Tt4 

'response = MsgBox(Msg7, Style, Title) 

Cells(22, 4). Select 

Msg7 = "The Tt4 lower limit (as indicated) is too low to converge, increase Tt4 lower limit and try again!" 
response = MsgBox(Msg7, Style, Title) 

'warn = 1 
GoTo 10 
End If 

If itbinit = 0 Then 
Tt4p5 = Tt4 * tau tH 

Call enthalpy(l, Tt4p5, 0, 0, Cpitb, gamaitb) 
tau alp itb = Cpitb * Tt4p5 / (Cpc * T o) '(11) 
fitb = 0# 
tau_itb =1# 

GoTo 5 
End If 

Call enthalpy(l, Tt4p5, 0, 0, Cpitb, gamaitb) 

25 

tau alp itb = Cpitb * Tt4p5 / (Cpc * T o) 

If tau_alp_itb <= (tau alp b * tau tH) Then 
itbinit = 0 

Tt4p5 = Tt4 * tau tH 

Call enthalpy(l, Tt4p5, 0, 0, Cpitb, gamaitb) 
fitb = 0# 
tau_itb =1# 

tau_alp_itb = Cpitb * Tt4p5 / (Cpc * T o) 

GoTo 5 
End If 

If itbinit = 1 Then 
’add MSH on 8/1/04 for ITB -1 
’Component: INTERSTAGE-TURBINE BURNER 
tempfitb = eta itb * hPRitb / h_0 - tau_alp_itb 
If temp fitb - 0# Then 
response = MsgBox(Msgll, Style, Title) 

’Msg7 = "Occurs at " & "PIc=" & pi_cH * picL & ", PIf=" & pi_f & ", Mo=" & M0 & ", alpha=" & alpha 
& ", Tt4=" & Tt4 

'response = MsgBox(Msg7, Style, Title) 

'If singlept <> 1 Then jump = MsgBox(MsgJump, Stylel, Title) 
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'warn = 1 
GoTo 10 
End If 

'Fuel model 1: CSH 

If Sheets("Input").ComboBox4. Value = "CSH" Then 
fitb = (taualpitb - tau_alp_b * tau_tH) / (tempfitb) * (1# + fb) '9/3/04 

'(22a) 

Else 

'Fuel model 2: MSH 
loop4 - 0 

Tt4p4 = Tt4 * tautH ’(22b) 

Call enthalpy(0, Tt4p4, fb, h_t4p4, 0, 0) 

f_4p5i = (1# + fb) * (tau alp itb - tau alp b * tau tH) / (temp fitb) 

'(22c), initial value for f_4p5i 


4 


If loop4 > 20 Then 
loop4 = 21 
GoTo 10 
End If 

Call enthalpy(0, Tt4p5, f_4p5i, h_t4p5, 0, 0) 
temp fitb = eta itb * hPRitb - h_t4p5 
If temp fitb = 0# Then 
response = MsgBox(Msgl 1, Style, Title) 

'Msg7 = "Occurs at " & "PIc=" & pi cH * pi cL & ", PIf=" & pi f & ", Mo=" & M0 & ", alpha=" 
alpha & ", Tt4=" & Tt4 

'response = MsgBox(Msg7, Style, Title) 

'warn = 1 

'If singlept <> 1 Then jump = MsgBox(MsgJump, Style 1, Title) 

GoTo 10 
End If 

fitb = (1# + fb) * (h_t4p5 - h_t4p4) / (temp fitb) '(22d) 

If Abs(fb + fitb - f_4p5i) > 0.0001 Then 
f_4p5i = fb + fitb 
loop4 = loop4 + 1 
GoTo 4 
End If 
End If 

'lkhl - added cp=cp(T,f) and gama=gama(T,f) on 1 1/15/04 
Errfitb = Abs(fitbi - fitb) 

If Err fitb > 0.0001 Then 
fitbi = fitb 

Call enthalpy(l, Tt4p5, fitb, 0, Cpitb, gamaitb) 

GoTo 25 
End If 
Tkh2 

'add MSH on 8/1/04 for ITB -2 

tau_itb = Tt4p5 / (T o * taur * tau c * tau_b * tau tH) '(24) 

End If 


5 

ritb = (gamaitb - 1#) * Cpitb / gamaitb '(3) 

f = fb + fitb '(23) 

tautL = 1# - (tau cL - 1# + alpha * (tau f - 1#)) * tau r / ((1# + fb + fitb) * tau alp itb * effm tL) 


'(26) 

pi tH = tau tH A (gamat / ((gamat - 1#) * e_tH)) '(27) 

pi tL = tau tL A (gamaitb / ((gamaitb - 1#) * e_tL)) '(28) 
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'(31) 

tempM9 = Pt9_P9 A ((gamaitb - 1#) / gamaitb) 

If tempM9 < 1# Then 

response = MsgBox(Msg8, Style, Title) 

'Msg8 = ’’Occurs at ” & ”PIc=” & pi_cH * picL & ”, PIf=” & pi_f & ”, Mo=” & MO & ”, alpha=” & alpha 
& ”, Tt4=” & Tt4 

'response = MsgBox(Msg8, Style, Title) 

'warn = 1 

'If singlept <> 1 Then jump = MsgBox(Msgsheet, Stylel, Title) 

GoTo 10 
End If 

Tt9_T0 = taur * tau_c * tau b * tau tH * tau_itb * tau tL '(32) 

T9 T0 = Tt9_T0 / Pt9_P9 A ((gamaitb - 1#) / gamaitb) '(33) 

M9 = ((tempM9 - 1#) * 2# / (gamaitb - 1#)) A 0.5 '(34)'5/5/04 

M8 = M9 

If M9 >= 1 Then M8 = 1 

v9_a0 = M9 * (gamaitb * ritb * T9 T0 / (gamac * re)) A 0.5 '(35) 
v9 - v9_a0 * aO 
v9_v0 = v9 / vO 

Ptl9_P19 = P0P19 * pi_r * pi_d * pi_f * pifn '(36) 

tempM19 = Ptl9_P19 A ((gamac - 1#) / gamac) 

If gamac < 1# Or tempM19 <1# Then 
response = MsgBox(Msg9, Style, Title) 

'Msg8 = "Occurs at " & "PIc=" & pi_cH * picL & ”, PIf=" & pi_f & ”, Mo=” & MO & ”, alpha=” & alpha 
& ”, Tt4=" & Tt4 

'response = MsgBox(Msg9, Style, Title) 

'warn = 1 

'If singlept <> 1 then jump = MsgBox(Msgsheet, Stylel, Title) 

GoTo 10 
End If 

M19 = ((tempM19 - 1#) * 2# / (gamac - 1#)) A 0.5 '(37) 

M18 = M19 

If M19 >= 1 Then Ml 8 = 1 

Ttl9_T0 - tau r * tau f '(38) 

T19 T0 = Ttl9_T0 / Ptl9_P19 A ((gamac - 1#) / gamac) '(39) 
vl9_a0 = M19 * T19 T0 A 0.5 '(40) 

vl9 - vl9_a0 * aO 
vl9_v0 = vl9 / vO 

fe me = (1# + f) * v9_a0 - M0 + (1# + f) * (ritb / re) * (T9 T0 / v9_a0) * (1# - P0 P9) * (1# / gamac) 

’(41) 

ffmf = vl9_a0 - MO + (T19T0 / vl9_a0) * (1# - P0P19) / gamac 

'(42) 

F_mo = (aO / gc) * (fe me + alpha * ff mf) / (1# + alpha) '(43) 
unitchange =1# 

If Sheets('Tnput”).ComboBoxl. Value = "SI” Then 
unitchange = 1000000# 

Else 
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unitchange = 3600# 

End If 

s = f/((l# + alpha) * Fmo) * unitchange 


'(44) 


'(44) 

If Fmo < 0 Or s < 0 Then wdata - 0 
'Efficiencies 

corenume = (1# + f) * v9 * v9 - vO * vO 'corrected on 6/16/04 


fanume = alpha * (vl9 * vl9 - vO * vO) 

nume = 0.5 / gc * (corenume + fanume) '(45) 

denu = fb * hPRb * eta b + fitb * hPRitb * eta itb '(46) 

etathermal = nume / denu * 100# '(47) 

eta_P = (1# + alpha) * F mo * vO / nume * 100# '(48) 

eta_0 = (1# + alpha) * F mo * vO / denu * 100# '(49) 


'Mach = 0.3 at HP turbine exit 
M4p4 = 0.3 

Tt4p4 - tau_r * tau_c * tau_b * tau_tH * T o 

T4p4 = Tt4p4 / (1# + 0.5 * (gamat - 1#) * M4p4 * M4p4) 

v4p4 = M4p4 * (gamat * rt * gc * T4p4) A 0.5 

dTitb = Tt4p4 * (tau itb - 1#) 

GoToll 

10 

'If warn = 1 Then 'added on 1 1/30/04 

Msg = "Occurs at " & "PIc=" & pi_cH * pi cF & ", PIf=" & pi_f & ", Mo=" & MO & ", alpha=" & alpha & 
Tt4=" & Tt4 

response = MsgBox(Msg, Style, Title) 

If last = 1 Then 
jump = 1 
GoTo 11 
End If 

If singlept <> 1 Then jump = MsgBox(Msgsheet, Style 1, Title) 

'End If 
11 

End Sub 
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Appendix B — Flowchart of Iterative Solution Scheme in 
Microsoft® Excel Code 


Main Code Structure 


ITB On-design 
Excel code 


Multiple-point 

Calculation 


Single-point 

Calculation 


Engine Station 
Data 


1 


2 


3 


4 


5 


Compressor 
Pressure 
Ratio (tt c ) 


Flight Mach 
Number (M 0 ) 


Fan Pressure 
Ratio (7t f ) 


Fan Bypass 
Ratio (a) 


Main Burner Exit 
Temperature (T t4 ) 


Figure B.l. — Flow-chart of the main structure for ITB on-design Microsoft® Excel code. 
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On-design Mlultiple-point Calculation 



Figure B.2. — Flow-chart of the iterative solution scheme for on-design multiple-point calculation. 


Remarks : *DLV - Do Loop Variables; **FLV - For Loop Variables. 
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On-design Single Point calculation 


Begin ^ 


Read user’s 
option 


z 


Read input data 


1 

' 

ITB computation 

1 

l 


Write input data 


c e "‘ ) 


7 


Figure B.3. — Flow-chart of the solution scheme for ITB on-design single point calculation. 
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Engine Station Data 


Beg inj) 


Z 


Read user’s option 


Read input data j 


Z 


1 

r 


ITB computation 


1 

F 


Read M 0 at each station 


7 


Compute pressure and 
temperature at each station 


z 


Write output data 


CZD 


Figure B.4. — Flow-chart of the solution scheme for 1TB engine station data computation. 
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Figure B.5. — Flow-chart of the iterative solution scheme for ITB computation subroutine. 
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Appendix C — Combustion Fuel Burn Models 

C.l Types of Fuel Burn Models 

Two models (p. 116 of ref. 4) are being used to calculate the fuel burned in the main burner and the 
ITB: 

(1) Constant Specific Heat ( CSH) model - The air and combustion gases at inlet and exit of each 
component are modeled as “calorically perfect gases” with constant specific heats. The values of the 
specific heats are different at inlet and exit of two combustors (main burner and ITB ) 

(2) Modified Specific Heat {MSH) model - All engines properties are calculated using CSH model 
except the fuel used. In this model, the inlet total temperature ( T t3 and T t4A ) of combustors are calculated 
using CSH model while the exit total temperatures ( T t4 and T tA5 ) are directly obtained from the user 
inputs. Nevertheless, the total enthalpies (h t 3 , h t4 , h t4A , and h tA5 ) will be calculated directly from the 
Variable Specific Heat ( VSH) model (p. 116 of Ref. 4). Provided h PR is given, the amount of fuel burned 
is then calculated using eq. (B4) and (B7), respectively. Therefore, the improvement on enthalpy 
calculation gives better estimates of fuel used. 

The following ENTHALPY subroutine is based on pages 105-106 of reference 3 and used in the MSH 
model to calculate the enthalpy (h) at engine stations t3, t4, t4.4, and t4.5. 

Subroutine Enthalpy 

ENTHALPY ( 0,T,f h,0,0) 

Inputs'. T and / 

Output : h 


h • = 

n air 


h ref +A 0 T + Et2 + — T 2 + — r 4 +^-T 5 +Llj 6 +^-T 1 +L-T% 1 (C.l) 

7 2 3 4 5 6 7 8 


^ prod 


href +A 0 T + ^-T 2 + —T 2 +^-T 4 +^-T 5 +^-T 6 +^-T 7 +^-T 8 1 (C.2) 

' 2 3 4 5 6 7 8 \ prod 


h • = 

n mix 


hair fh prod 

TT7 ' 


(C.3) 


where the constants A 0 , etc. are given in Table C.l (p.106 of Ref. 3). 


Table C.l.— Constants for air and combustion products used in subroutine ENTHALPY(ft T, f h, c v ,y). 


Constant 

Air alone 

Combustion products of air and (CfDn fuels 

^0 

2.5020051 x 10~‘ 

7.3816638 x 10 _1 

A i 

-5.1536879 x 10 -5 

1.2258630 x 10 -5 

A 2 

6.5519486 x 10~ 8 

-1.3771902 x 10~ 8 

A 3 

-6.7178376 x 10“ 12 

9.9686793 x 10' 12 

A 4 

-1.5128259 x 10“ 14 

-4.2051 104 x 10“ 14 

A 5 

7.6215767 x 10 -18 

1.0212913 x 10~ 18 

a 6 

-1.4526770 x 10“ 21 

-1.3335668 x 10“ 21 

a 7 

1.0115540 x 10“ 25 

7.2678710 x 10“ 25 

href 

-1.7558886 Btu/lbm 

30.58153 Btu/lbm 
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The equations needed to calculate the fuel-air ratios of the main burner (f b ) and the ITB (f, b ) are 
listed below for both the CSH and MSH models. 


Main burner 
If CSH model then 

p ^X-b Z^cLScH_ 

r\bhpR / h 0 - Tx~b 

Else 

7] 3 = Tq ^r^cL^cH 

ENTHALPY (0, T t3 , 0, h t3 , 0, 0) 

Set initial value of fuel/air ratio at station 4 =f 4 \ 

3 ENTHALPY (0, T t4 , f 4i , h t4 , 0, 0) 

p _ ~ ^3 

r\b hpR ~ h , 4 


(C.4) 


(C.5) 


(C. 6 ) 


If \fb-Ui\> 0.0001, then / = fa and go to 3; else continue. 
End if 

Inter-stage Turbine Burner 
If CSH model then 


fab =(i+4/> V SCf" 

l\itb "PR / «0 t'k-itb 


Else 


Tt 4.4 - T t MtH 


ENTHALPY (0, T t4A , f b , h t4 4 , 0, 0) 

Set initial value of fuel/air ratio at station 4.4 =fa.s { 
4 ENTHALPY (0, T t4 . 5 , f 4 _ 5i , h t 45 , 0, 0) 


fitb -(l + fb) 


^4.5 ~^?4.4 
T1 itbhpR ~h t 4.5 


If | fb + fitb ~ / 4 . 5 /I > 0.0001, then fa 5i = fa + fab and go to 4; else continue. 


End if 


If fab 0 then y it b — Yt , Cp^b — Cp ^ , and ^ 


(C.7) 


(C. 8 ) 
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